Quantum and thermal fluctuations in the SU(N) Heisenberg spin-glass model near the 

quantum critical point 
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We solve for the SU(N) Heisenberg spin-glass in the limit of large N focusing on small S and T. We 
study the effect of quantum and thermal fluctuations in the frequency dependent response function 
and observed interesting transfers of spectral weight. We compute the T— dependence of the order 
parameter and the specific heat and find an unusual behavior for the latter at low temperatures 
in the spin-glass phase. We find a remarkable qualitative agreement with various experiments on 
the quantum frustrated magnet SrCr9pGai2-9pOi9. 



PACS Numbers: 75.50.Lk, 75.40. Gb, 75.10.Jm 

Disordered quantum magnets are fascinating systems. 
The understanding of the interplay between disorder, 
quantum and thermal fluctuations remains among the 
most challenging problems of condensed matter physics 
[1-4]. These three aspects are always present to some 
extent in experiments on real systems, therefore a clear 
understanding of their interplay is very desirable. In 
systems where disorder is relevant we usually encounter 
the phenomenology of slow dynamics that is associated 
with glassy states. When quantum fluctuations become 
important the phases with glassy orders can be driven 
to more conventional phases through interesting quan- 
tum phase transitions [5]. One example that is captur- 
ing the interest of experimentalist and theorist alike is 
LiHo2;Yi_a;F4 which is a dipolar coupled random mag- 
net [6] and has been recently the focus of beautiful ex- 
periments [7] where quantum fluctuations are introduced 
and controlled by means of a transverse magnetic fleld. 
An other example, and perhaps the archetype of frus- 
trated quantum magnets, is the bi-layer Kagome lattice 
SrCr9pGai2-9pOi9 (SCGO) that only becomes a spin- 
glass at the low temperature of about 5K. This compound 
has been thoroughly investigated over the years [8-12] 
and, in sharp contrast to ordinary classical spin-glass sys- 
tem, exhibits some unusual remarkable features that are 
associated with strong quantum fluctuations: The mag- 
netic fluctuation spectrum, x"{^)i is found to vanish lin- 
eraly in oj at low frequencies [10] and the speciflc heat 
is proportional to [9]. On the theoretical side, these 
observations have remained largely unaccounted for. 

The progress in the understanding models of disor- 
dered quantum magnets in finite dimensions is rather 
slow. In fact, a great deal of our knowledge still re- 
lies on solutions of systems with long-ranged interac- 
tions. These mean-field models are appealing because 
they are mathematically more tractable while retaining 
much of the physics associated with slow dynamics. It 
is worth pointing out that in many actual systems, such 
as LiHo2;Yi_a;F4 that is an insulator, the magnetic in- 
teractions do have power-law decay, thus each individ- 



ual spin interacts with others well beyond their nearest 
neighbours [6]. 

Among the simplest mean field models for quantum 
spin-glasses, the quantum version of the Sherrington- 
Kirkpatrick (SK) model received a great deal of atten- 
tion. It is a Heisenberg model with gaussianly distributed 
random interactions between all pair of spins in the lat- 
tice. The model was first considered by Bray and Moore 
[13] and they predicted a spin-glass phase at low tempera- 
ture, substantially reduced from the usual (Ising) version 
of the SK model. Further progress was prevented be- 
cause replica symmetry broken solutions were expected 
at low T. Later, Sachdev and Ye introduced a generaliza- 
tion of the model to SU(N) spins which could be studied 
in the large N limit [14]. They found a very interest- 
ing spin liquid phase down to T = 0. In more recent 
work on this model, a generalized phase diagram as a 
function of T and S was obtained using a bosonic repre- 
sentation [15-17]. The spin quantum number S can be 
thought of a parameter that controls the strength of the 
quantum fluctuations. For S— s- oo one goes to the "clas- 
sical" limit while for small S the quantum fluctuations 
are strongest. A low temperature spin-glass phase was 
found for all non zero S and Tg ~ at large S [15-17]. 
Remarkably, the spin-liquid phase was also found at very 
low spin S [14-16]. Therefore, quantum fluctuations can 
drive the model through an interesting quantum critical 
point between a spin-liquid state at S— >0 and a quantum 
spin-glass for finite S. 

Recent numerical studies based on quantum Monte 
Carlo [18] and exact diagonahzation [19,20] techniques 
for the SU(2) model have validated some aspects of pre- 
vious investigations. 

The goal of the present work is to focus on the different 
roles played by quantum and thermal fluctuations in the 
SU(N) SK model within the quantum critical regime. We 
obtain the detailed behavior of the dynamical spin sus- 
ceptibility for small S and T, both in the paramagnetic 
(PM) and spin-glass (SG) phases. We flnd interesting 
transfers of spectral weights in the magnetic response. 
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We also find that the spin-glass order parameter has a 
simple temperature behavior at small S and obtain the 
correct specific heat at low temperatures. In addition, 
wc discuss the remarkable qualitative agreement that we 
find between our model solutions and the experimental 
results in the SCGO compound that we mentioned above. 
The model Hamiltonian is 
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(1) 



where the magnetic exchange coiiplings Jy arc indepen- 
dent, quenched random variables distributed according 
to a Gaussian distribution where J is the variance and 
the unit of energy. As already pointed out by Bray and 
Moore [13], one uses the replica trick to average over 
the disorder [2] and the lattice infinite-range model maps 
exactly onto a self- consistent single site model with the 
action (in imaginary time r, with /? the inverse temper- 
ature) : 



Seff = Sb- 



2N 



r drdr' Q''\t - t')'^\t) ■ ^'(r') 
Jo 



(2) 



and the self-consistency condition 



Q'^\t - r') = ^ < ~§\r) ■ ^\r') >s^„ (3) 



where a,b = 1, • • • , n denote the replica indices (the limit 
n ^ has to be taken later) and Sb is the Berry phase 
of the spin [14]. Due to the time-dependence, the solu- 
tion of these moan-field equations remains a very difficult 
problem for N = 2, even in the paramagnetic phase [18]. 

We shall use the bosonic representation [14-17] for the 
spin operators where S is represented with Schwinger 
bosons b by Sai3 — b\bp — SSaf), with the constraint 
^^6t^6„ ^ SN {fd < S). In the language of Young 
tableaux, these representations are described by one line 
of length SN. They are a natural generalization of an 
SU{2) spin of size S. 

In the N ^ oo limit, the mean field self- 
consistent model (2-3) reduces to an integral equa- 
tion for the Green's function of the boson G'^''{t) = 
— < T6°(T)6t''(0) > where the bar denotes the average 
over disorder and the brackets the thermal average [14] : 



Sf(T)= J^Gt\T))'Gt\-T) 



(4) 



(5) 



Gr(r = 0") = -S 

The local spin susceptibility XIoc{t) 
is given in the large N limit by XIoc{t) = 



(6) 

=< S{t)S{0) > 
Gl%T)Gr{-T) 




FIG. 1. The imaginary part of the dynamical spin suscep- 
tibility x"('-^) as a function of w for various values of S at 
T = 0.04. 



In the spin glass phase it is enough to perform a one 
step symmetry broken solution [15-17]. Equations (4-6) 
were solved self-consistently on the Matsubara axis. To 
obtain the imaginary part of the w-dependent dynami- 
cal response x"oc('^)' ^li'- solutions were analytically con- 
tinued to the real axis using a method based on Pade 
approximants [21]. The general form of the spin suscep- 
tibility can written as x'Ld'^) = ^^i^^{'^)+x'reg{'^) where 
qea is the spin glass order parameter. 

In Fig. 1 we show results for x"eg('^) ^■t low T=0.04 
and several values of S across the PM-SG boundary. At 
this T, the critical S is found at S «0.28. We observe a 
qualitative change in the regular part of the response as 
S is increased. At low S, in the PM phase, the suscepti- 
bility shows the finite temperature spin liquid behavior, 
obeying x"('^) ^ tanh{ui/2T) for u) small. On the other 
hand, as S increases and the system goes into the SG 
phase and x'regi'^) opens a pseudogap. The thermal ex- 
citations become gradually less important and a linear in 
uj behavior shows up clearly. In fact, we find that the low 
frequency behavior is proportional to u/S (inset). This 
is consistent with the large S solution obtained in [15]. 

Wc now turn to the role of thermal fluctuations for 
fixed S = 1/2 that enables comparisons to numerical 
results obtained in the SU(2) model [18-20]. The freez- 
°^ ing temperature is found at Tg ssO.133 in good agree- 



ment with all previous estimates [13,18-20]. We start 
at low T = 0.05 well in the spin-glass phase. In Fig. 2 
we show the susceptibility with clean pseudogap ~ u) be- 
havior at low CO (the 5{u)) is not shown for clarity). As 
T is increased, one observes that excitations gradually 
fill the pseudogap with a narrow low frequency feature 
that peaks at w ^ 0(T). These excitations come from 
the gradual melting of the frozen spins, i.e., from the 
^-function part. Another interesting effect that one ob- 
serves is that spectral weight from high frequencies of 
order J is transferred down to the pseudogap. The in- 
terpretation of this is that when spins are frozen in the 
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FIG. 2. The imaginary part of the dynamical spin suscep- 
tibility x"{<^) as a function of uj for various values of T at 
S = 0.5. Below Tg~0.133 the systems is in the SG phase 
and x"{^) develops a delta function part at a; = that is not 
shown for clarity. 

spin glass state, they still have a fast motion of preces- 
sion around the axis of their local frozen field. That 
motion originates a contribution to the susceptibility at 
Lu ^ 0{J) with a strength proportional to the frozen frac- 
tion, i.e., to the order parameter q^;^. As T increases, 
the spins (and thus their local field) melt, so the contri- 
bution from the motion of precession gradually decreases 
and merges with the excitations of order T that now fill 
the pseudogap. We may point out that this behavior is 
qualitatively similar to the results obtained from exact 
diagonalization of small SU(2) clusters [20]. 

As T is further increased one enters the PM phase and 
the melted peak fully merges with the higher frequency 
part of x"i^) ^iid there is no more a clear separation of 
energy scales. In this quantum disordered regime the low 
frequency behaviour of x"(<^) is oc tanh{u)/2T) as in the 
spin-liquid state [14]. 

Thus we have seen that the regular part of the response 
begins at low T in the SG phase with a clean and linear in 
Lo pseudo-gap, then above Tg the; gap becomes thermally 
filled down to very low frequencies, and finally when T is 
well above Tg the pseudo-gap clears up again displaying 
once more a linear in lo behavior. We find remarkable 
that this unusual evolution is qualitatively identical to 
that reported from neutron experiments in SCGO (cf. 
Fig. 3 of Ref. [10]). Moreover, the neutrons have also 
revealed that the spatial correlations are extremely short 
('^ 2.5A) which may render additional justification to the 
relevance of the present mean-field theory results. 

In order to better appreciate the evolution of the trans- 
fers of spectral weight at low frequencies and for small 
T, it is useful to consider the spectral density, defined 
by p{uj) = x"oc('^)/(s~'^'^ ~ 1) t^^t obeys the sum-rule 
/ Preg{i^)du + qea = S{S + 1). The spectral density at 
different temperatures is shown in Fig. 3 where the inten- 
sity of the delta function part is denoted by the height of 




FIG. 3. The regular part of spectral density function p{u)) 
as a function of w for various values of T at S = 1/2. Below 
Tg ~0.133 the systems is in the SG phase. The inset shows 
p{u)) in the PM phase (left), entering the SG phase (center) 
and at low T within the SG (right). The height of the arrow 
indicates the spectral intensity of the delta part at the origin. 

the arrow (see inset) and can be thought as the fraction 
of frozen spins. In the main panel we show the regular 
part of the spectral density that has a broad background 
contribution that remains almost temperature indepen- 
dent. Most of the T— dependence occurs at the low fre- 
quencies where a rather narrow peak is present. At the 
higher temperature, in the PM phase there is no delta 
function contribution, however the large peak at small u) 
indicates that a portion of the degrees of freedom actually 
got slowed down (left panel of inset). As T is lowered, the 
system enters the SG phase and the peak becomes nar- 
rower and losses weight. A (5-function contribution thus 
emerges as some of the slow spins become frozen (central 
panel of inset). When T is further lowered towards T = 
we observe how the resonance losses all its weight that 
gets transferred to the 6 part (right panel of inset). The 
strong quantum fluctuations are responsible for the large 
remanent background spectral density that corresponds 
to a large fraction of spins remaining disordered. 

The results for the spin-glass order parameter qEA{T) 
arc shown in the inset of Fig. 4. The behavior at moder- 
ate and large values of S was previously investigated in 
Rcf. [16]. We focus here in the small S and T regime and 
find that the order parameters obeys the simple form 

qEA{T) = qEA{T = 0) - aT^ (7) 

with a a constant, which is similar to the solution of the 
SK model [1,22]. We also find that the value of qea at 
T = and its jump at Tg{S) are strongly reduced by 
quantum fluctuations when S ^ 0. In fact, the later 
vanishes faster than in contrast to the quadratic de- 
pendence at large S [16]. It is interesting to note the good 
agreement between qEA{T = 0) « 0.20 for 5* = 1/2 with 
the corresponding estimate qEA{T = 0) w 0.18 obtained 
in an exact diagonalization study of the SU(2) model [20] . 
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FIG. 4. The specific heat Cv{T)/T for S = 1/2 and 1. 
Inset: The SG order parameter [qEA{T = 0) - qEA{T)]/T as 
a function of T for the same values of S. 

Possibly more significant is to mention that a frozen mo- 
ment of only 20% of the total possible elastic value was 
observed in neutron scattering on SCGO [10,11]. 

Finally we turn to the specific heat Cy . In earlier work 
[15], it was argued that the specific heat is linear in T 
at low temperatures within the SG phase. We calculated 
the specific heat in the whole range of S values, with 
particular care in the small S regime. In contrast to 
previous results [23], we find that the behavior within 
the SG phase at low T has the form Cy{T) oc as 
is shown in Fig. 4. This result can be justified by the 
following argument: The energy of the system is given 

by, 

E{T) = J^^ dTQ-\T)Q^\-T). (8) 

Writting this expression in the SG phase in terms of the 
one-step replica symmetry broken solution [15] and in- 
serting the functional form of the order parameter (7), 
one finds that the lowest order terms in the T— expansion 
of E{T) are proportional to T"^, and T^. It is not 
hard to check that the coefficient of the divergent 
term exactly cancels. On the other hand, the linear term 
should also vanish as otherwise the groundstate would 
have extensive remanent entropy. Therefore, one may 
expect that the first non vanishing contribution to the 
energy at low T is proportional to T^, which implies a 
quadratic behavior for C„ (T) as is in fact borne out of our 
calculations. While ordinary spin-glasses usually show a 
linear in T specific heat, we find remarkable that the un- 
usual behavior that we obtain was indeed observed 
among the early investigations of the spin-glass state of 
the SCGO compound [9] that pointed to the important 
role that quantum fluctuations might play in that com- 
pound. 

To conclude, we have solved for the behavior of the 
Heisenberg SU(N) spin- glass in the limit of large N and 



investigated in detail the parameter region of small S and 
T close to the quantum critical point of the model. We 
observed the qualitatively different role played by quan- 
txvai and thermal fluctuations through their effect on the 
transfers of spectral weight in the dynamical response 
functions. We obtained the functional form of the or- 
der parameter and present an argument to support the 
flnding of an unusual quadratic behavior of the speciflc 
heat at low temperature within the spin glass phase. We 
flnd a very interesting qualitative agreement with vari- 
ous experimental findings in the SCGO compound one 
of the most thoroughly investigated quantum spin-glass 
system. Extensions of our work to incorporate a more 
realistic geometric structure might be interesting routes 
for future research. 
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